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Abstract. Weyl's unitary operators for displacement in position and momen- 
tum commute with one another if the product of the elementary displacements 
equals Planck's constant. Then, their common eigenstates constitute the Zak 
basis, with each state specified by two phase parameters. Accordingly, the trans- 
formation function from the position basis to the Zak basis maps the Hilbert space 
on the line onto the Hilbert space on the torus. This mapping is one-to-one pro- 
vided that the Zak basis states are periodic functions of their phase parameters, 
but then the mapping cannot be continuous on the whole torus. 
With the periodicity of the Zak basis enforced, the basis has a double Fourier 
series. The Fourier coefficients identify a discrete basis which complements the 
periodic Zak basis to form a pair of mutually unbiased bases. The discrete basis 
states are the common eigenstates of the two complementary partners to the two 
unitary displacement operators. These partner operators are of angular momen- 
tum type, with integer eigenvalues, and generate the fundamental rotations of the 
torus. Conversely, the displacement operators are the ladder operators for their 
partners. 

For each consistent phase convention for the periodic Zak basis, and thus for 
the line-onto-torus mapping, there is a corresponding discrete Zak basis and a 
corresponding pair of partner operators. Examples of particular interest are the 
conventions that give a continuous mapping in one phase parameter or are sym- 
metric in both phase parameters. The latter emphasizes the Heisenberg— Weyl 
symmetry between position and momentum. 

We discuss briefly the relation between the Zak bases and Aharonov's modular 
operators. Finally, as an application of the Zak operators for the torus, we men- 
tion how they can be used to associate with the single degree of freedom of the 
line a pair of genuine qubits that are potentially entangled. 
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1. Introduction 

A convenient choice of variables may render simple a complex physical problem. 
Joshua Zak introduced [1,2] what he termed "the kq representation" as a particularly 
convenient mode of description for electrons in solids [3,4]. Zak's work was based on 
Weyl's [5] presentation of unitary displacement operators for coordinates and momenta 
(obtained by exponentiating both the momentum and the coordinate operators) and 
noting that for some particular choice of parameters these displacement operators 
commute. Weyl's work [5] led to several studies that dealt with both mathematical 
and physical representations wherein the spatial coordinates arc defined modulo some 
conveniently selected length a and the momenta modulo 27r/a (in units for which 
h = 1). However, it was Zak [1,2] who was the first to recognize the wide context of 
the approach and to study it in a systematic way [6]. 

The Zak transform of a state may be regarded as a mixed position-momentum 
(or time-frequency) representation of the state. The transform enabled Zak [3, 4] to 
reduce the problem of an electron in a periodic potential in a constant magnetic field 
to a form that made possible studies of further physical applications — to digital data 
transmission, for instance, as is considered in [7]. 

Reference [6] is a thorough review of the Zak transform and its applications up 
to the year 1988; for recent applications, see [8-10] for example. In the more recent 
literature the representation finds its widest use in signal processing where time plays 
the role of the spatial coordinate and frequency that of the momentum; see [7, 11, 12] 
and references therein. 

In the present work, which we regard as a contribution to our understanding of 
quantum kinematics, we consider the Zak transform as a mapping of the Hilbert space 
on a line to the Hilbert space on a torus. This allows us to study, in sectional periodic 
and discrete Zak bases which are mutually unbiased [13,14]. In section|21we use some 
freedom in the permissible definition of the phase of the periodic wave functions in 
the Zak basis to present and discuss three convenient phase choices. Section 0] is 
devoted to the study of the discrete Zak basis; we display and discuss the Wigner 
functions [15] (adopting the conventions of [16]; see the footnote at l|4.5|) below) that 
are associated with the wave functions in these bases. In section [S] the position and 
momentum operators are given as differential operators in the Zak representation of 
section [3J Section deals with the relation of the Zak bases and operators to the 
modular variables introduced by Aharonov et al. in Ref. [17]. A possible application 
for the Zak bases in quantum information theory is hinted at in section We close 
with a Summary. 

2. Setting the stage 

We consider a single continuous quantum degree of freedom, the motion along the 
x axis, for which the position operator X and the momentum operator P constitute 
the fundamental complementary pair of hermitian dynamical variables. They obey 
the Heisenberg-Born commutation relation [X, P] = ih, and we adopt the usual 
normalization conventions for their complete sets of eigenbras and eigenkets. 

We denote Weyl's unitary operator for the momentum shift P — > P — po by U 
and that for the position shift X — > X — xq by V, 

U = e~ ipoX / n , v = e ixoP / h . (2.1) 
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When both are acting on an operator, a function of X and P, their order is irrelevant, 
f(X, P) - f(X -x ,P-po) = UWf(X, P)UV 

= V^U\f(X,P)VU . (2.2) 

But the order does matter for the transformation of bras and kets, as illustrated by 
(x\ UV = e-^/^ (x\ V = c -'^ x / n (x + x \ , 

(x\VU= (x + x \U = e - i K>(* + *<>)/» {x + x \ , (2.3) 

unless the product pqXq is an integer multiple of Planck's constant h = 2irTi. Following 
Zak [1,2], we therefore choose the two displacements such that 

p x = 2irh, UV = VU. (2.4) 

The product poxo then equals the phase space area that is roughly associated with 
one quantum state, and it is fitting to regard po and xq as corresponding atomic units 
for momentum and position. This is emphasized, for example, by the relations 

U = e -2vriX/a;o ^ y _ c 2niP/p ^ 

2ni[P/p 0l X/x ]=l. (2.5) 

Now that U and V commute, they have common eigenstates, which we label by 
the real phases a and (3 that are associated with the unit-modulus eigenvalues, 

U\a,0) = \a,P)e~ ia , 

V\a,(3) = \a,/3)e i P. (2.6) 

Since the characterizing eigenvalues are 27r-periodic in a and /?, we require that the 
eigenstates are periodic as well, 

\a, P) = \a + 2tt, (i) = \a,0 + 2tt) . (2.7) 

In the completeness relations, then, 

d/3 



I 6a f 

J(2tt) 27r J(2tt) 27r 



\a,0)(a,0\=l, (2.8) 



the integration is over any 2ir interval, and the wave functions of the Zak 
representation, ip(a,f3) = (a,(3\ ), are functions on the torus, not on the square. 
In the orthonormality relation 

(a,f3\a',(3') = {2^f 5 {2 ^ {a - a')S (27I \f3 - /?') (2.9) 

we meet the periodic delta function 

S (2 *\v)= E %-27rfc) = — elkiP - ( 2 - 10 ) 

h— — ■ oo k— — oo 

It appears also in the explicit construction of the projector on an a, (3 state in terms 
of the unitary shift operators, 



oo - . 

\a,(3)(a,p\ = ]T (e ia Uy (e'^v) 



(2ir) 2 5 (27T \p Q X/h - a)5 {2 ^{xvP/h - (i) . (2.11) 
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The Fourier coefficient kets \l, m) that are defined by the double Fourier series of 
\a,(3), 



|a,/3>= £ \l,m)e- [ V a - m P\ (2.12) 

l.rn— — oc 

form the discrete Zak basis that is normalized in accordance with 

oo 

(Z,m|Z',m') = 5 mm >Su> , ^ \hm){l,m\ = 1. (2-13) 

Lrn — — oo 

Upon observing that 

(a,/3|/,m) =e i ( /a - m ^, |(a, /3|Z, m)\ = 1 , (2.14) 

we note that the periodic and the discrete Zak bases are mutually unbiased. 

We regard the discrete basis states \l, m) as common eigenstates of two integer 
operators, 

L\l,m) = \l,m)l, M\l,m) = \l,m)m. (2.15) 
Their action on the continuous Zak states amounts to differentiation, 

\ d c) 
( a ,/?|L= T — < a ,/3| , (a,p\M = i— (a,/3| , (2.16) 
i da. op 

which implies that they are the generators of the fundamental rotations of the torus, 

(a, p\ e i( - a ' L ~ = (a + a',/3 + f3'\ . (2.17) 

Conversely, the unitary displacement operators are ladder operators for the discrete 
Zak basis, 

U l> V m ' \l, m) = \l- I', m - m') . (2.18) 

The commutation relations 

[U,V]=0, [U,L]=U, [U,M]=0, 

[V,M]=V, [M,L]=0, [V,L]=Q (2.19) 

reveal that the pairs (U, L) and (V, M) are the dynamical variables of two independent 
degrees of freedom of azimuth-angular-momentum type. These are, of course, the two 
degrees of freedom of the torus. 

3. Wave functions: Periodic Zak basis 

As a consequence of the eigenket equations (|2.6|l . the position wave function {x\a, 0) 
of the periodic Zak state, which is the transformation function between the position 
representation and the periodic Zak representation, obeys the functional equations 

(x\a,0) =S a -P° x ^(x\a,l3) 

= c-' ll3 {x + x a \a,(3) . (3.1) 

The corresponding equations for the momentum wave function are 

(p\a,/3) =e ia (p + P oM) 



c 



- xop/h) 



(p\a,(3). (3.2) 
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Their general solution is 



(x\a,/3) = ^x(a,(3)e- i ^e i ^ X0 S ( - 2 ^{p x/h^a), 



(p\a,0) =^ x (a,P)e i ^e- [a P/P°5^(xop/h-(3), (3.3) 



where x( a > P) must be of unit modulus, 

\x(a,P)\=l, (3.4) 

for consistency with the normalization (|2.8I) , and the periodicity of (|2.7[) requires that 
X(a, (3) obeys 

X (a, (3) = X (a + 2tt, p)e~^l 2 = e ia / 2 X (a, (3 + 2tt) . (3.5) 

Any X { a i0) permitted by these constraints is an acceptable phase convention for the 
periodic Zak bases. 

Equation (|3.5(l implies that 

X (a, 13) = (-l) ab e i( - al3 ~ ab ^ 2 X (a - 2™, (3 - 2nb) (3.6) 

if a, b are integers, so that x( a > P) is specified by stating its values within a standard 
2-7T x 2ir square. We find the following three choices for x(a,/3) particularly natural 
and interesting: 

(a) X {a,(3) =e iQ (^" L&1) , 

(b) x(a,f3) = e _i/3 fe~ L^l) , (3.7) 

(c) x(a, /9) = (-1) L^l L^l e 1 ( L^l a " t L^l ) , 

where 

|/|= / dC <5(C-fc-i) = z arctan(tan(7rz)) (3.8) 

A:— — oo 

denotes the integer closest to z. The x( a j/^) of choice (a) is continuous in a but 
discontinuous in /?; that of choice (b) is discontinuous in a and continuous in (3\ 
and the one resulting from choice (c) is symmetric in a and (3 in the sense of 
%(a, (3) = x{~f3, ol) = x(/3, a)* and discontinuous in both. 

As these examples demonstrate, discontinuities are a generic feature of the wave 
functions we are dealing with here, and of other functions derived from them. It 
will be necessary on various occasions to assign values to such functions at their 
discontinuities. We employ the convention to assign half the sum of the one-sided 
limits, for which 

!+l for x > 0, 
for x = 0, 
-1 for x < 0, 

[sgn(a;)] 2 = l for aU x, (3.9) 
e'V for x > 0, 

^sgn(x) = } CQSV fora . = , 

e~ i{ P for x < 0, 
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are illustrative examples. As an application, and for future reference, we note that 



1 d 



(a) 



47T 



2tt 



X(«, PT r ^-X(«, /3) = <! (b) - £ + <5 (2,r) - tt) sin , 



i 9a 



Atx 



(c) - 



2tt 



+ 5 (27r) (a-7r)(-l)L^l sin f 



'(a) - +5 (27r) (/3-7r)sina, 

47T 



(3.10) 



x(a,/3)*i^gX(a,/3) = < 



(b)-- 

1 
2 



(c) 



Q 

2^ 
a 
2^ 



+ ^ 27r )(/3-7r)(-l)L^l 



sm 



2 ' 



for the respective choices in 1(3.7(1 . 

Further, as an important consequence of (|3.5|l . we have the two Fourier series 



■ a/3 . 

x(ot,0)=e \ A(a — 2irm)e 



im/3 



m=—oo 
oo 



= e 4tt \ 



(3.11) 



l= — oo 



with 



(27T) 



— e 4* x(a,/3) , 

Z7T 



A(a) 

J(2ti-) 



(3.12) 



where the integrands are 27r periodic functions of the respective integration variable. 
Either one of the single-argument functions A(a) or /i(/3) introduced here, which are 
related to each other by Fourier transformation, 



d(3 i2l 



*(«)=/ ^e^/?), M(/3)=/ — e" 1 ^A(a), (3.13) 



2tt 



da 



2vr 



can thus be used for a complete specification of \( a jP)- 

The restriction imposed by the normalization 13.4(1 can be stated as 

da 



2vr 



e" Q A(a - 2vrm) e 1( ' a A(a - 27rm') 



(3.14) 



or equivalently as 



oo 

[e ito A(a - 2ttto) e" Q 'A(a' - 2nm) * = 2ir5(a - a') (3.15) 



L,m= — OD 



or by the corresponding equations for Equation l|3.14(l has the appearance of 

an orthonormality relation, while ((3.15(1 looks like a completeness relation, which is, 
of course, not accidental and will be clarified in section 0] 
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For the particular choices of l|3.7jl . we have 
(a) \(a) 



2 . a . fa 
— sin — = sine — 
a 2 \2 



1 for — 7T < (3 < 7T 

else 



360 : 



(b) A(a) = S a Q , fj,([3) = sine 



/3 



(3.16) 



(c) \(a) = sine 



a + 2na 



where a, b are the integers a 



a 
2^ 



/3 + 2nb 



and 6 = 



2- 



We emphasize the symmetry of choice (c), where A (7) = /i(7) is its own Fourier 
transform. Figure^shows a plot of this function for — 87r < 7 < 87r. 



1 



0.8 



0.6 



0.4 



0.2 



-0.2 




-8 -6 



6 8 



Figure 1. The discontinuous function A(7) = fi(-f) of choice (c) in equation 
13.161 which is its own Fourier transform. 



4. Wave functions: Discrete Zak basis 

The wave functions for the discrete Zak basis states \l, m) are obtained from the wave 
functions in (|3.3JI by Fourier analysis. The position wave function 

<x|l,m>=/ pf f( X \a^)e^-^) 

J(2tt) ^ J(2tt) 111 

1 e UP°*/ h \(p oX /h - 27rm) (4.1) 
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reveals the physical significance of X(a) of equation l(3.11jl . and we recognize now 
that (|3.14|) and (|3.15|1 are the position representations of the orthogonality and 
completeness relations in ((2.13|) . 

The corresponding momentum wave function 

(p\l m ) = J— e - imx ° p / n ti(x p/n - 2tt0 (4.2) 

involves /i(/3), so that the Fourier relations of (|3.13l) are actually those between the 
position and momentum representations. We note the consistency with Ij2.18|l . which 
is quite explicitly visible in 

\l,m) = U- l V- m \0,0) : (x\l,m) = e ilpoX / n (x - mi„ |0, 0) , 

(p\l, m) = e - imx °P/ n (p - lp \0, 0) , (4.3) 

with 

(x|0,0) = X( Po x/h)/^, (p|0,0) - n(*oP/K)/y/Po- ( 4 - 4 ) 
For a visualization of the states of the discrete Zak basis, it may be helpful to 
plot their Wigner functions [15]f, given by 

Wi m (x, p) = Woo (x - mx Q ,p-lpo), 

W 00 (x,p) = 2 J ^ X( Po x/h + a)* e i2a P/P° X( Po x/h - a) 

= 2 J^L n(x oP /h - /?)* <PM** »(x p/h + f3) , (4.5) 

so that Wi m is centered at (x,p) — (mxojpo) if Woo is centered at (x,p) — (0,0), as 
is the case for the three choices in (13.711 or (13.16(1 . Choice (b) yields 

f 2snc(l - 2\x/x \,2np/p ) for |x| < x /2 , 
W$\x,p) = { (4.6) 
{ for |x| > x /2, 

where 

snc (z, a) = Sm(za) = - [ dC e <a for < z < 1 (4.7) 

is the incomplete sine function, sinc(a) = snc(l, a). For choice (a) the roles of x and 
p are interchanged, i.e., 

W$ (x, p) = (x oP /p 0lPo x/x ) . (4.8) 

Figure □ shows W$(x,p) for -0.75 < x/x < 0.75 and -5.25 < p/p < 5.25. 
For choice (c) we obtain 

W$(x,p) = 2 (-l) {a + Mb + k) snc (\l-\t\-k\, (2a + s+j S gn(s))^) 



i,fc=0,l 



|l-| s |-j|,(26 + t + fesgn(t))^) , (4.9) 



where a = [2x/xq] and b = l2p/po\ are the integers closest to 2x/xq and 2p/po, 
respectively, and s,t account for the differences, 

2— = a + s, 2— = b + t with - - <s,t<-. (4.10) 
x Po 2 2 

I dix dp 

f They are normalized in accordance with / W(x,p) = 1 and bounded by —2 < W(x,p) < 2. 

/ 2nh 
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-0.5 -0.3 -0.1 0.1 0.3 0.5 

x/x 



Figure 2. The Wigner function Wq \x,p) corresponding to the ket |0,0) for 

choice (b) in equation 13.161 . Top: above the x,p plane, with the quadrant 

for x > 0,p < cut out. Bottom: The WjJq = contour lines where the Wigner 

function changes sign (solid lines), and the Wqq^ = 1 contour line where the 
Wigner function is at half-maximum value (dashed line) . 
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-2-1012 

X/X 



Figure 3. The Wigner function Wqq (x,p) corresponding to the ket |0, 0) for 
choice (c) in equation 13.161 . Top: Wqq above the x,p plane, with the quadrant 

(c) 

for x > 0,p < cut out. Bottom: The = contour lines where the 

Wigner function changes sign (solid lines), and the Wqq = 1 contour line where 
the Wigner function is at half-maximum value (dashed line). The dotted square 
identifies — 1 < x/xq,p/pq < i the central peak region of = 0. 
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fc) 

The interchange pqx <-> xop has no effect on Wg , 

W&W) = W$(x p/p ,p x/x ) , (4.11) 
which is yet another manifestation of the symmetry that characterizes the phase 
convention (c) of ((3.7(1 . 

In effect, then, the phase space is tiled by squares of size \xopo = ^nh, centered 
at {x,p) = (\axo, \bpo), with the pair s,t parameterizing the points of each square 
relative to the center. At the centers of the tiles we have 

2 if a = and b = , 



W$ (±axo, \bp ) = 2<5 a0 <5, 



60 



ifa^0or&^0. 



(4.12) 



Figure |3 shows Wqq\x,p) for —2.25 < x/xq,p/pq < 2.25, the vicinity of the maximum 
at (x,p) = (0,0). 

These plots demonstrate that the Wigner functions in 1(4.6(1 and ((4.9(1 are 
continuous in both x and p, although their marginals — of which 



dpW(x,p) = p \\(p x/h)\ 
dxW(x,p) — xo\(jl(xop/Ji)\ 



(4.13) 



are special cases — exhibit discontinuities as a rule. Indeed, it is a matter of inspection 
to verify that and Wqq are continuous. 

5. Operator relations 



We combine the wave functions of 14.1(1 and 1(4. 2JI with the Fourier series of 1(3.11(1 to 
establish 
(x\ e iaL ~ & M \p) 



(x\p) 



y^ e -ix P /n Y^{x\l, m)S la ~ m ®{l, m\p) 



fp x x p\ 



exp i 



. axop — Ppqx 



fp x x p \ 



The normalized matrix elements of L and M axe then obtained by differentiation, 



{x\L\p} p 1 d 

\x\p) 2p i da 

( x \ M \p) x d 

\x\p) 2x dp 



a = pox/h ' 
f} = xop/h 

a = pox/h ■ 
j3 = x p/h 



(5.2) 



The required derivatives are available in ((3.1U|) . and so we arrive at 



(a)L = 



P 
Po 
P 



M 



•TO 



V h 



sm ■ 



(b) L = S ( ' 

Po 

1 I P 



1 fx 

M = -{- 



Po 
X 

x 



iw) ( P^_ 

\ n 
p 

Po 
X 

Xq 



— 7r sm ■ 



x P 



M = 



PoJ£ 

h 

x' 

Xq 



(5.3) 



_ 5 (2n) ( XOP 



sm ■ 



x P 

2h 
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which express L and M in terms of X and P for the three phase conventions of (|3.7|1 . 

The reverse relations, which state the line variables X and P in terms of the 
torus variables U, L, and V, M, are at hand as soon as one observes that L — P/po and 
M — X/xo are 2ir periodic functions oipoX/h and xqP/U, that is: they are functions 
of U and V. For the three choices of (|3.7|l we have 

'(a) S {21T) (P-Tr)sma, 



X 



= M+ < 



a 

1 / a 



(b) 



a 
2^ 



2^ 



,,1-a.l . a 
l)l2n I sm — : 
2 



= L 



(a) 



2tt 



2tt 



(5.4) 



(b) (5( 27r >(a-7r)sin/3, 

) +<5( 2Tr )(a-7r)(-l)L^l sin | 





/3" 


) 


\2tt 


_2tt 





where e ia = {/ and e'' 3 = are understood in the doubly periodic functions on the 
right, as illustrated by 



, „ . -s-1 . a .fa ai\ v-^ 

(-1) I sm - = sm - tt [— )= ^ 

(-1)* 



sin(ja) 



j = -oo 

For the example of choice (a), then, we have 

1 a 



^ i(2j + 1)tt v ; 



(5.5) 



(or,/3| X = x ( -y^g ( a >/?l + sin " ("^ 



2tt 



(5.6) 



i da 2-7T 

for the differential-operator representation of X and P in the continuous Zak 
representation. 



6. Aharonov's modular operators 

Adopted to the present notational conventions, Aharonov's modular position and 
momentum operators [17], X m and P m , are introduced by the symmetric relations 

X = x N x +X m , P = Po N p + P m , (6.1) 

where N x and N p are the integer operators 

N x = [X/xo] , N p = LP/pol • (6.2) 

In view of 
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we can regard the modular operators as the logarithms of the unitary operators U and 
V in the sense of 



ih e e + U 
X m = X - xo [X/xq] = — In 

Po e e + U 1 

h e e + V 
Pm = P-Po[P/Po] =— In- 



> e -+ 



(6.4) 



ip e s + V - > e _> 

where the limiting procedure ensures — ^xq < X m < \xq and — \po < P m < \po- 

Therefore, (X m . N p ) is a pair of complementary observables that is equivalent 
to (U,L); and, likewise, the pair (P m ,N x ) is as good as (V, M). But the two pairs 
(X m , N p ) and (P m , N x ) together do not refer to the two independent periodic variables 
of the torus because N p does not commute with N x , 

i[N X: N p ] =l--S^\p Q X/h~TT)-S^\x P/h-TT). (6.5) 

This non-commutativity is as it should be, because the pair (X m , N p ) belongs to the 
phase convention (a) with N p — L^ a \ whereas the pair (P m ,N x ) goes with phase 
convention (b), as N x = M (b >. 

The non-commutativity of N x and N p reminds us that it is not possible to 
constrain a quantum system both to a finite x range and to a finite p range. It 
appears that, in the present context of the Zak bases, the symmetric compromise of 
imperfect simultaneous localization is given by the discrete Zak basis states of the 
symmetric phase convention (c). Indeed, Wqq (x,p) is very strongly peaked within 
the central square — ^^x/xo,p/po^^ that is indicated in Figure 



7. Hinting at an application: Toroidal qubits 

The three operators that are introduced by 

tn+io-a = [l + {-l) L ]U, <7 3 = (-l) L (7.1) 

are such that 

tr{cri} = tr{cr 2 }=tr{o-3} = G, 

2 2 2 1 

°1 = a 2 = CT 3 = 1 . 

eri<T2=io" 3 , 0-2C3 = ifi , ay,ai—\a 2 , (7.2) 
so that they are the components of a genuine Pauli vector operator a = {a\, 02, 03). 
We can, therefore, use them to associate a qubit with the U, L degree of freedom. 
Likewise, 

ri+ir 2 = [l + (-l) M ]F, r 3 = (-l) M (7.3) 

defines a second toroidal qubit for the V, M degree of freedom. 

For any quantum state on the line, there is then a two-qubit state with the 
statistical operator 

p=-(l + a-(a) + (f) - T + a- (af) ■ r) . (7.4) 

The two vectorial expectation values (ct), (t) together with the dyadic expectation 
value (ctt) constitute the 15 parameters that specify the two-qubit state. Once their 
values are known, one can, for instance, decide whether p is separable or not, and 
thus whether an entangled qubit pair is concealed in the given state on the line. This 
application and others are, however, beyond the scope of the present paper on the 
properties of the Zak bases. 
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8. Summary 

The unitary operators for displacement in position and momentum commute with one 
another when the product of the elementary displacements forms the elementary area 
h = 2-kU in phase space. The Zak basis is composed of the complete set of eigenstates 
of these commuting operators and is therefore specified by two phase parameters. 

Hence the Zak basis maps the Hilbert space on the line onto the Hilbcrt space on 
a torus. We showed that for this mapping to be one to one, periodicity of the Zak basis 
states in their phase parameters is required and this, in turn, imposed discontinuity in 
the phases. We used the freedom present for the definition of this phases discontinuity 
to consider three choices which we deemed most natural: (i) one phase parameter 
continuous, say a, and the other, j3, discontinuous; (ii) conversely, a discontinuous 
while (3 continuous; and (iii) symmetric discontinuity in both phases. We thereby 
generalized Zak's approach, whose original choice [1-4] was asymmetric. 

The resultant periodic Zak basis may be represented by a double Fourier series. 
The Fourier coefficients constitute a doubly discrete representation basis that is 
mutually unbiased with the Zak basis. The states of this discrete representation are 
characterized by pairs of integer eigenvalues. The hermitian operators, for which these 
states are eigenstates with these eigenvalues, generate the fundamental rotations of 
the torus. Thus with each consistent phase mapping convention of the line onto the 
torus we have one pair of mutually unbiased Zak bases. 

Further, we considered briefly the relation between the Zak bases and Aharonov's 
modular operators and interpreted the operators employed by the latter in terms of 
the former: both, per force, reflect the basic fact that it is not possible to constrain a 
quantum system to finite domains in both position and momentum. 

As a possible application of our study of the Zak operators on the torus we 
mentioned the possibility of associating with a single degree of freedom on the line a 
pair of potentially entangled qubits. This is achieved by a rather simple, and quite 
natural, construction of two sets of Pauli operators in terms of the basic unitary 
operators of the Zak bases. 
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